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Abstract 

In diffractive jet production, two high energy hadrons A and B collide and 
produce high transverse momentum jets, while hadron A is diffractively scat- 
tered. Ingelman and Schlein predicted this phenomenon. In their model, 
part of the longitudinal momentum transferred from hadron A is delivered 
to the jet system, part is lost. Lossless diffractive jet production, in which 
all of this longitudinal momentum is delivered to the jet system, has been 
discussed by Collins, Frankfurt, and Strikman. We study the structure of 
lossless diffractive jet production in a simple model. The model suggests that 
the phenomenon can be probed experimentally at HERA, with A being a pro- 
ton and B being a bremsstrahlung photon with virtuality Q 2 . Lossless events 
should be present for small Q 2 , but not for Q 2 larger than 1/Rp, where Rp 
is a characteristic size of the pomeron. 
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I. INTRODUCTION 



In high energy pp collisions at hadron colliders, one sometimes produces high transverse 
momentum jets in the final state. Here high transverse momentum means, say, transverse 
momentum greater than 20 GeV. These jets are interpreted as the decay products of quarks 
and gluons (partons) produced by the hard scattering of two partons from the incoming 
hadrons. Cross sections for 



calculated in QCD according to this picture match the observations well [0,0. Normally the 
two hadrons A and B are broken up in this process. It is clear that this should be so, since 
a 20 GeV momentum transfer is quite sufficient to disrupt a hadron. Nevertheless, in 1985 
Ingelmann and Schlein || predicted that events of the type 



where hadron A is diffractively scattered, should occur with a small but not tiny probability. 
Here by "diffractively scattered," we mean that A emerges with a fraction (1 — z) > 0.9 of 
its original longitudinal momentum and with a small transverse momentum |P^| < 1 GeV. 
The transverse momentum transfer can also be characterized using the invariant momentum 
transfer t from the hadron: t = (P A - Pa) 2 = -(P' A 2 + z 2 M\)/(l - z) « -P' A 2 . 

The picture for such diffractive hard scattering proposed by Ingelman and Schlein is 
that hadron A exchanges a pomeron with the rest of the system, where "pomeron" means 
whatever is exchanged in elastic scattering at large s, small t. Thus the cross section is 
proportional to the pomeron coupling to hadron A as measured in elastic scattering. The 
pomeron carries transverse momentum — and a fraction z of the hadron's longitudinal 
momentum. Here we do not need to know what a pomeron is, only that its momentum 
is carried by quarks and gluons. One of these collides with a parton from hadron B to 
produce the jets. Let the parton that participates in the hard scattering carry a fraction 
x of the longitudinal momentum of the incoming hadron A, and thus a fraction x/z of 
the longitudinal momentum transferred by the pomeron. Then the cross section in this 
model is proportional to a function f a /p(x/z,t; /i), where f a /p(£,,t;fi)d£, is interpreted as 
the probability to find a parton of kind a in a pomeron, where the parton carries a fraction 
£ of the pomeron's longitudinal momentum. In principle, f a /p can depend on the invariant 
momentum transfer t, and it should also depend on a scale parameter /i that characterizes 
the scale of virtualities or transverse momenta in the hard parton scattering. 

Several guesses for the £ dependence of / a /p(£,t;/i) were given in the literature || ||. 
but ultimately it was left to experiment to measure these functions. An analysis of the 
theoretical expectations and of the approximate factorization hypothesis that is inherent in 
the model of Ingelman and Schlein was given in || . 

The reaction (|2|) anticipated by Ingelman and Schlein has been seen at the CERN collider 
by the UA8 experiment ||. However, the experiment suggests a feature not anticipated in 
0-0. It was expected that the functions f a /p(x/z) would have support only for x < z. That 
is, some of the momentum fraction x transferred from hadron A would be lost, appearing 
in low Pt particles instead of the jets. Instead, the experiment suggests that a fraction of 



A + B -> jets + X 



(1) 



A + B ->• A + jets + X, 



(2) 
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the events are lossless in the sense that x = z. It is as if the formula for the cross section 
contained a term proportional to 6(1 — x/z). 

It will be important to confirm the existence of a term that is proportional to 6(1 — x/z) 
or, more generally, that is singular as x/z — > 1. The UA8 experiment demonstrates that 
something is present in addition to terms that vanish as x/z — > 1. However, given the 
experiment's rather limited resolution, it is not clear to us that the experimental findings 
could not be fit with a function that is merely finite as x/z — > 1. We shall present evidence 
in a forthcoming paper that such terms are to be expected [|7| . 

Recently, the Zeus and HI experiments at HERA have reported the first evidence for the 
analogous reaction in deeply inelastic electron scattering ||, 

e + A -> e + A + X. (3) 

Here the diffractively scattered proton A has not been seen, but there is evidence for the 
"rapidity gap" expected in such events. Information on the x and z dependence of the cross 
section is not yet available. We will discuss HERA physics later in this paper, concentrating 
one + A— > e + A + X events in which the virtual photon has low virtuality, rather than 
high virtuality. 

The purpose of this paper is to consider some of the theoretical issues in light of the 
experimental results. In particular, we reexamine the factorization hypothesis of 0, and 
argue that this hypothesis is at best approximate. This conclusion is in agreement with that 
of Collins, Frankfurt, and Strikman ||, although our analysis of the source of the breakdown 
of Ingelman-Schlein factorization will be presented in more depth than that of Ref. ||. We 
frame much of the analysis in terms of a simple model that allows a detailed examination 
of the theoretical issues. 

Our analysis suggests three conclusions. First, the factorization model with a "distribu- 
tion of partons in a pomeron" may ultimately not be a fruitful way of understanding the 
process. Second, the lossless events found in hadron-hadron collisions should be absent in 
deeply inelastic scattering. 

Third, perhaps the most interesting conclusion suggested by the model is that if one 
examines the process e + A^e + A + Xasa function of the virtuality Q 2 of the virtual 
photon emitted by the electron, then the lossless events should be present for Q 2 = but 
disappear gradually as Q 2 becomes larger than some value 1/Rp. The length scale Rp may 
be interpreted in the model as the transverse "size" of the pomeron. 

Our discussion of the experimental possibilities at HERA bears some resemblance to that 



found in a recent paper of Brodsky, Frankfurt, Gunion, Mueller, and Strikman JT0[. These 
authors study e + A — >e + A + V + X, where A is a proton and V is a vector meson. They 
look at the region 1 GeV 2 < Q 2 C s. As in the present paper, large Q 2 implies a small 
size for the wave function of the quarks produced by the virtual photon. As in the present 
paper, these quarks couple to the proton via two gluons 



II. DIFFRACTIVE JET PRODUCTION 

Let us consider, as a definite example, the inclusive cross section for the production of 
two jets in a high energy collisions of two hadrons, A and B. Let the initial hadron A have 
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momentum 

M 2 



P£ = (P+, P A , P A ) = (P+, ^r, 0) , (4) 



where we denote P ± = (P° ± P 3 )/v2 and where we denote transverse vectors by boldface 
letters. Similarly, hadron B enters the scattering with momentum 

PS = (Pb + , Pb, Pb) = (^L, Pb" 0) • (5) 

We specify the two jets by variables Ex, Xa, and Xb, given in terms of the four momenta 
Pi and P<2 of jets 1 and 2 by 

P T = (|P 1 | + |P 2 |), 
X A = (P 1 + + P 2 + )/P A + , 

x fl = (pr + p 2 -)/p s - . (6) 

The production of high transverse momentum jets is a hard process, to which QCD 
factorization applies. Thus the total cross section for inclusive two jet production can be 
written as 

da(A + B -> jets + X) 
dEx dX A dXs 

f , . , s f , . , \ da(a + b — > jets + X) 
2^ y d %a fa/A{x a ; V) J dx b f b/B (x b ;n) dX ' O 



a,b 



Corrections to this asymptotic equality are suppressed by a power of m/Ex, where m rep- 
resents a typical hadronic momentum scale, say 300 MeV. Here da is the parton-level cross 
section, computed in perturbation theory with a certain prescription for removing infrared 
divergences. The functions f a /A(x a ',fi) are the parton distribution functions. Notice that 
the Born-level cross section da is proportional to 5(x a — Xa) 5(x b — Xb) since, at the Born 
level, all of the momentum carried by the colliding partons is transferred to the observed 
jets. In higher orders of perturbation theory, one will have Xa < x a and X B < x b . 

The topic of this paper is inclusive jet production in which we make a very restrictive de- 
mand on the final state. We demand that hadron A appear in the final state with momentum 
P' A given by 

The invariant momentum transfer to hadron A is 

t = (P A - Pa) 2 = ~yz~ z { P a + z " m2 } ■ ( 9 ) 

For diffractive scattering, we require that t be small compared to s = (Pa + Pb) 2 and that 
z be small, say z < 0.1. Since Xa < x a < z, one must look for diffractive hard scattering in 
the region X A < 0.1. We thus consider the cross section 

da diS (A + B ^A+ jets + X) 



dEx dXA dXs dz dt 

in the region described. 



(10) 
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III. MODEL OF INGELMAN AND SCHLEIN 



We begin our discussion of diffractive jet production by reviewing the model of Ingelman 
and Schlein ||. We follow much of the notation of Ref. ||, but break the presentation 
into two stages. In the first stage, we hypothesize that the cross section for diffractive jet 
production can be written in terms of a diffractive parton distribution: 

da diS (A + B -> A + jets + X) 

dEx dXA dXg dz dt 
^ f , d fa/A( x a, fi) f , , . . da(a + b -> jets + X) 

g J dXa dzdt J dXb M X « A dE T dX A dX B ■ (U) 

Here 

dzdt dXa {U) 

represents the probability to find in hadron A a parton of type a carrying momentum fraction 
x a , while leaving hadron A intact except for the momentum transfer (z,t). In the second 

diff / 
%/A\ 



stage, we hypothesize that d /^(i a , /x) / dz dt has a particular form: 



Here j3(t) is the pomeron coupling to hadron A and a(t) is the pomeron trajectory. The 
function f a / P (x a /z,t, fi) is then the "distribution of partons in the pomeron." 

In writing Eq. (0), one thinks of the pomeron as a continuation in the angular momen- 
tum plane of a set of hadron states. Since hadrons contain partons, the pomeron should 
also. Thus one has in Eq. (|T3|) the standard factors describing the coupling of the pomeron 
to hadron A, together with a distribution of partons in the pomeron Eq. ([TT|) , on the 

other hand, is more basic. It says only that factorization still applies when hadron A is 



diffr actively scattered. It is Eq. (11) that will be of concern in this paper. We shall argue 



that, in fact, factorization does not generally apply when hadron A is diffractively scattered. 



IV. FACTORIZING AND NON-FACTORIZING GRAPHS 

In this section, we discuss the Ingelman-Schlein model in terms of the contributing 
Feynman graphs. We first give a graph-based argument that the Ingelman-Schlein picture 
is plausible. Fig. [I] shows a typical graph that contributes to diffractive jet production. The 
shaded circles represent Bethe-Salpeter wave functions for the two hadrons, A and B. In the 
center of the diagram, there is a tree level contribution to the hard scattering of two gluons 
with momenta approximately collinear to the momenta of their respective parent hadrons. 
The hard scattering produces the observed high Et jets. The jets carry plus-momentum 
X^Pa- The gluon that enters the hard scattering from hadron A carries plus-momentum 
XaPjt, which equals X^P^ f° r Born-level hard scattering. We may view this gluon as a 
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normal part of the gluon cloud that accompanies any hadron. Since only a small fraction 
Xa < 0.1 of hadron A's plus-momentum has been lost to the hard scattering, it seems 
plausible that the gluon cloud and the constituent quarks can reassemble themselves into 
the hadron bound state. It is only necessary to dispose of one more gluon, as indicated in 
the figure. Then the net color transfer from hadron A can be zero. If this gluon carries a 
small momentum fraction x' into the final state, then the net momentum fraction lost from 
hadron A, z — Xa + x', will also be small (z < 0.1, say). Similarly, the net transverse 
momentum — transferred from the hadron can be small (|P^| ~ 1 GeV). 

diffjet.epsf here: just uncomment the macro. 

FIG. 1. Diffractive jet production. 

Let us now compare Fig. [I] with a typical graph that contributes to pomeron exchange, 
i. e. hadron-hadron elastic scattering at large s and small t, as depicted in Fig. 0. Of course, 
the pomeron is much more complicated than this, so the reader is invited to imagine some 
much more complicated graphs Nevertheless, it seems an attractive proposition that 
whatever graphs contribute to the bottom half of the pomeron in Fig. ^| can contribute 
equally to the bottom half of Fig. [D. 

pomeron. epsf here: just uncomment the macro. 

FIG. 2. Typical graph contributing to the pomeron. 



This graphical argument makes the Ingelman-Schlein model seem plausible. Since the 
pomeron-proton coupling is large, the argument also suggests that the cross section for 
diffractive jet production at a given is not an infinitesimal fraction of the total cross 
section for jet production with the same jet Et- In this paper, however, we will focus on 
two other features of Fig. [I], factorization and longitudinal momentum loss. 

First, we note that the top half of the graph is connected to the bottom half only at the 
hard scattering. The absolute square of the graph breaks up into a convolution of a factor 
that contributes to the distribution function of gluons in hadron B, times a contribution to 
the hard scattering cross section, times a contribution to the distribution function of gluons 
in hadron A, with the additional requirement that hadron A is diffractively scattered. This 
is just the structure of the factorized cross section formula ([11]). Thus the usual factorization 
program (as described in |T2J) is relatively straightforward for this graph (although there 
are some subtleties related to the gluon polarizations). 

Second, we note that, in order to get the color right, it was necessary to emit a gluon (or 
more than one gluon) into the final state. Inevitably, then, some of the net plus-momentum 
zP + transferred by the pomeron is not contributed to the hard scattering, but is lost into 
the final state in the form of low transverse momentum particles. Thus the contribution to 
the cross section from graphs like Fig. [l] will be nonzero only for Xa/ z < 1. 

Consider, however, the graph of Fig. |3|. Here the gluon that had to be emitted in order 
to make the color right is absorbed by a parton from hadron B, which is moving in the 
opposite direction. As we will see in the following sections, the plus-momentum transferred 
by the extra gluon in this process is negligible compared to zP + . Thus the contributions 
from graphs like Fig. |3| will be effectively proportional to 5(1 — X^/z). Furthermore, graphs 
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like this break the Ingelman-Schlein factorization represented by Eq. (|TT|). As argued in 
|T2|1 , such factorization breaking effects cancel when one calculates an inclusive cross section 



to produce jets with no restriction on the final state. However, here there is an important 
restriction on the final state: we allow only final states consisting of an elastically scattered 
proton. Thus factorization breaking effects may survive. As noted in 0, a signature of 
these effects is the appearance of 5(1 — Xa/z) terms in the cross section. 

diffjetmod.epsf here: just uncomment the macro. 

FIG. 3. Diffractive jet production with color transferred to spectator parton. 

In the following sections, we explain a simple model in which this effect can be examined 
in some detail. Within this model, we find that the factorization breaking effect has some 
simple features. We abstract the simple features from the model and propose that one can 
test for these features in electron-proton scattering at HERA. 



V. A SIMPLE MODEL 

We will examine diffractive jet production within the context of scalar-quark QCD. In 
this model, there is an SU(3) gauge field A£(x) as in normal QCD. There is also a color 
triplet quark field qi(x) with quark mass m, but we take qi(x) to be a scalar field instead of 
a Dirac field. Since the theory includes a scalar quark field, a 4-quark coupling is necessary, 
but we set the renormalized coupling constant g± to a negligibly small value. The lagrangian 
is 

C = (D„q)\D»q) - mV? - \g%G^ - ±(^)(^) 

+ Faddeev-Popov terms — ^(q^q) 2 - (14) 

This theory has the same behavior as spinor QCD for collinear and soft gluon emission from 
quarks. Its chief advantage is that it allows a perturbative model for a quark- ant iquark 
bound state. We introduce a scalar, color singlet meson field <f>(x) and couple it to the 
quarks using 

C^ = G<P(x)q\(x)q t (x). (15) 

We work to lowest nontrivial order in the (pq^q coupling G, letting the <pq^q vertex play the 
role that is played by the (amputated) Bethe-Salpeter wave function of a meson in spinor 
QCD. We denote the mass of the meson by M and take M to be smaller than 2m, so that 
the meson cannot decay into a quark and an antiquark. This model has been used for similar 
purposes in Refs. and |L4[ . 



In the following sections, we examine diffractive jet production within the context of low 
order perturbative graphs in this model. The idea is to understand some of the basic physics 
by means of examples. We will be working with perturbation theory for the cross section 
at order a^G 6 and a b s G & . Of course, the real strong coupling a s (fj) is not small except at 
large /i, and some parts of our graphs contain soft momentum flows, for which small \i would 
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be relevant. We hope (but do not prove) that the basic lessons we learn here will be valid 
at order a^G 6 for any N in the model, and at order in real QCD. However, such an 
analysis at arbitrary order N remains a challenge for the future. 



VI. WAVE FUNCTIONS 

The analysis that follows will make use of the null-plane wave function of the meson state 
in our model. For a meson moving in the minus direction, the wave function ip(x, k)$j is 
the amplitude to find that the meson with momentum P M = (M 2 /(2P~), P~ , 0) consists of 
a quark and an antiquark of colors i and anti-j respectively, with the quark having minus- 
momentum k~ = xP~ and transverse momentum k. The wave function is measured by 
operators defined on the null-plane y~ = 0. The precise definition, following the formalism 
of [|H^T7|, is 



k)« = 2x(l - x)P~ I d A y 5(y~) (0\ gi (y) qj(0)\P) . (16) 



Here we have chosen the normalization 

where Pi is the probability, which is of order G 2 , that the meson state consists of a {q,q') 
pair. In terms of the covariant (f>qq^ Green function amputated on the 0-leg, the definition 
(fL6|) can be written as 

r dk + 

^(x, k)« = 2x(l - x)P- / — G(k a , P% , (18) 
At lowest order in a s and G one has, 

Q(k a ,P% 3 =iGb l3 - - — — -. (19) 

k 2 — m 2 + it (P — k) 1 — m 2 + xe 

By integrating according to Eq. (|IBD , we find 

Gx(l-x) 
k 2 + m? - x(l - x)M 2 



= ^ I ^H V (2°) 



Notice that |f/>| 2 oc 1/k 4 for large k 2 . This good behavior in the ultraviolet, which arises 
from the fact that G has dimensions of mass, is the essential reason for the usefulness of this 
model. 

The wave function can be used to calculate, to order zero in a s , the probability for finding 
a quark in a meson i.e. the parton distribution function: 



3G 2 x(l-x) 
167r 2 m 2 — x(l — x)M 2 



|2 



Ki I 



(21) 
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We pause here to consider a more realistic example, chosen in anticipation of the applica- 
tion to HERA physics in the final section of this paper. We let the meson be a photon with 
polarization vector e and we let the quarks be genuine spin-1/2 quarks with negligible mass. 
The perturbative Green function for the photon to turn virtually into the quark- ant iquark 
system is 

ffM = fcCj ^L e . 7 ^^L. (22) 

Here e is the proton charge and Q is the quark charge in units of the proton charge. The 
Dirac indices are left implicit and there is an implicit unit matrix in the 3x3 color space. 
The wave function, defined analogously to Eq. flTB]) , is 

1 r fih+ 

tf(s, *) = ^j^ U xP~, k; P?) TV , (23) 

where U and V are the Dirac spinors for the quark and antiquark, respectively. 

We consider an off-shell photon with P 2 = —Q 2 < and zero transverse momentum, so 
that 

P> = {-■§=, P-,0). (24) 

To define the possible polarization states, we use light-like axial gauge with gauge-fixing 
vector f M = (1, 0, 0). In this gauge the photon propagator is iH^ j[P 2 + ie) where 

A/7- = _ r + + = Y: e(P, iYe{P if - e(P, L)^(P, L) v . (25) 

" ' V i=l 

Here e(P, l)' 4 and e(P, 2)^ are the two transverse polarization vectors, with v ■ e = P ■ e = 
and e • e = — 1. The longitudinal polarization vector is 

e(P,L)" = -Q-v». (26) 
P ■ v 

As Q 2 — > 0, only the transverse polarization contributes. The longitudinal mode becomes 
important for Q 2 ^ 0. It does not propagate, but is responsible for the "instantaneous scalar 
photon exchange" force [fT5| , |T6[] . 

We now insert polarization vectors into Q and perform the fc + -integration in Eq. (f2"B"D. 
For transverse polarization, we find 

^ k) = " 2P= + x(l - x)Q 2 ' (2?) 

where kj, = (0,0, k). This is quite similar to the scalar wave function, Eq. (EDI). The chief 
difference is the factor of k in the numerator, which leads to a logarithmic divergence in the 
normalization integral for if). (The spinors U and V depend on k, but this dependence is 
eliminated when the spinors stand next to a 7~.) For longitudinal polarization, we obtain 

^,k) = px(i-x)q k2 + ^ )Q2 . (28) 

Again, U'-f'V is independent of k. This wave function is small compared to that for trans- 
verse polarization when Q 2 <C k 2 but becomes comparable when Q 2 ~ k 2 . 
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VII. LOWEST ORDER CONTRIBUTION 



We consider the cross section 

da diS (A + B jets + X) 

dEx dXA dXs dz dt 



(29) 



for diffractive jet production in meson-meson collisions within the model described above. In 
Fig. |] we show a lowest order graph for this process. The jets are produced in a quark-gluon 
collision. The active quark in this collision comes from meson B. The gluon comes from 
meson A, where it has been emitted from one of the quark lines. Another gluon is emitted 
into the final state and meson A is reconstituted. There are eight ways of attaching the 
gluons to the scalar quark and antiquark lines; only one of them is shown. Similarly, there 
are four lowest order graphs for quark-gluon scattering. Thus one needs to sum over the 32 
Feynman graphs contributing to the amplitude. There is also an analogous set of graphs in 
which the hard collision is between an antiquark and a gluon. 

goodgraph.epsf here: just uncomment the macro. 

FIG. 4. Factorizable graph for diffractive jet production. 

The contribution to the cross section calculated from graphs like Fig. f| takes the factor- 
ized form of Eq. (|TTD in the high energy limit. This is a simple application of the general 
analysis of ||12|| . We do not give this analysis in any detail, but we do mention three key 
points. 

First, we discuss the polarization of the gluon carrying momentum k^ in Fig. |j. We 
consider all graphs to be in Feynman gauge in this paper, since to use a physical gauge 
introduces non-physical singularities into the graphs. In Feynman gauge, the relevant part 
of the graph has the form 

Ra Z^jP (30) 

where H represents the hard scattering part of the graph and J represents the meson wave 
function part at the bottom of the graph. Now we can replace — g a @ by 

E z a (kt, k4, A) k A , A) - " A + ^ + «V , (31) 



A=l 



where u a = (0,1,0) is a lightlike vector in the minus-direction and e a is the polarization 
vector for a gluon with momentum {k\, k^) and polarization A in null-plane gauge, u ■ A = 
A + = 0. (See ref. [I5[ for definitions.) This much is an identity. Now using the power 



counting of ref. [12|, one finds that the terms proportional to u a k A and u a u^ are negligible 
compared to the terms proportional to e a e f3 in the leading integration region, in which k A 
is approximately collinear to Pa- Thus these terms can be dropped. On the other hand, for 
each Feynman graph, the term proportional to k^u^ is larger than the terms proportional 
to e a e 13 . This is an artifact of Feynman gauge. When one sums Hak^ over the graphs 
contributing to H, there is an almost complete cancellation and the remaining terms are 
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negligibly small. The result of this is that — g a/3 can be replaced by J2 eC " ^ That is, only 
physical polarizations contribute. 

The second point concerns the gluon carrying momentum q in Fig. |j. The leading 
integration region for q is the collinear region, in which q is the same order of magnitude 
as the transverse momenta inside the meson: q 2 ~ m 2 . In an abelian theory, the graphs 
without this gluon would constitute the lowest order graphs for diffractive jet production. 
Fig. |] would represent another contribution, a contribution that is formally of higher order 
in a s but is still a leading contribution in the sense of not being suppressed by a power 
of ui/Et- However, in the non-abelian theory, Fig. ^ gives the lowest order contribution 
because the emission of a second gluon is necessary to restore the meson to a color singlet 
state. 

The third point concerns the factor representing meson B: 

Taking into account that P B = (M 2 / (2P£), P B , 0) and that (P B — k B ) 2 = m 2 , this factor is 

G6ij k 2 B + m 2 -x B (l-x B )M 2 ' (33) 



where x B = k B /P B and kg is the transverse part of k B . Comparing to the definition ( |20|) 
of the meson wave function ip, this is 

— lj){xB, k|) y . (34) 

xb 

Thus when we square the amplitude and integrate over kg we obtain the parton distribution 



function for meson B, Eq. (|21 ) 



In summary, the graph of Fig. |4] factors in a simple way. In the next section, we examine 
graphs that do not factor. 



VIII. COLOR EXCHANGE WITH SPECTATOR PARTON 

The graph in Fig. |5| is similar to that in Fig. |], but here the second gluon is absorbed 
on the spectator antiquark line from meson B rather than emitted into the final state. As 
in Fig. |4], there are eight ways of attaching the gluons to the scalar quark and antiquark 
lines associated with meson A; and there are four lowest order graphs for the hard quark- 
gluon scattering. Thus one needs to sum over the 32 Feynman graphs contributing to the 
amplitude. 

We will find two important differences between the cross section associated with Fig. |] 
and that associated with Fig. §. First, there is an important observable signature associated 
with Fig. Essentially all of the plus-momentum P£ — P'^ = zPj[ transferred from meson 
A appears in the system of high jets. A negligible amount is transferred to the spectator 
antiquark. Second, the cross section corresponding to graphs does not take the factored 
form flTTp. The cross section is not simply proportional to the probability to find the active 
partons from mesons A and B, but also contains information on the spectator interaction. 
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spectator. epsf here: just uncomment the macro. 

FIG. 5. Graph with spectator parton interaction. 



Consider the integration over in Fig. [5[ Using power counting arguments similar 
to those in Ref. jnj, one finds that the dominant integration region for the transverse 
components of g M is q 2 ~ m 2 . For q 2 much larger than this, the propagators in the loop for 
meson A are thrown far off shell, while for q 2 much smaller than this there is a cancellation 
between attachments to the quark and antiquark lines in meson A, since meson A is a color 
singlet. For q~ , the dominant integration region is q~ ~ m 2 /P^~, since, for q~ much larger 
than this the propagators in the loop for meson A are again thrown far off shell. Similarly, 
for q + , the dominant integration region is q + ~ m 2 / P B since, for q + much larger than this 
the propagators for the quark and antiquark in meson B are thrown far off shell. This 



integration region is denoted the Glauber region in Ref. [0 . 

We now examine more carefully the integration over q + . In the Glauber region, the 
only two propagators that have a significant dependence on q + are the spectator-antiquark 
propagator and the active-quark propagator. These have the form 

I ~ ! fQC) 

(P B -k B - q) 2 -m 2 + ie -2(1 -x B )P B q+ - A 2 + ie V ; 

and 

(36) 



(k B + q) 2 — m 2 + ie 2x B P B q + — A 2 , + ie ' 
where x B = k B / P B and 

A 2 = (k B + q) 2 -k|, 

1 - x B 



Notice that A 2 and A 2 are each of order m in the dominant integration region, in which x B 
is of order 1 and all k. B and q 2 are of order m 2 . 

The crucial fact about these propagators is that the poles in q + are on opposite sides of 
the integration contour along the real q + axis, creating a pinch singularity as discussed in 
T2]] . We can put the integral in a useful form by deforming the contour into the upper half 



q + plane, past the pole from the spectator propagator. The deformation can be continued 
until q + is large. Then the propagators are far off shell, so that the contribution from the 
deformed contour is negligible. The contribution that survives comes from the residue at the 
spectator pole. Thus we make the replacement q + —>■ — A 2 /[2(l — x B )P B ] in the active quark 
propagator. This value of q + is very small, and all other factors in the diagram depend only 
weakly on q + , so we may replace q + by everywhere else in the diagram. Thus we make 
the replacements 



[k B + q) 2 — m 2 + ie 



—i 



-A 2 + Al 



1 - x B 

-i (1 - x B ) 



(k B + q) 2 + m 2 - x B (l - x B )M 2 



(38) 
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followed by 

' 2n5(q + ). (39) 



(P B - k B - q) 2 - m 2 + le 2(1 -x B )P£ 
Finally, the interaction of the Glauber gluon with the spectator quark is 

-ig(2P B -2k B -q)H a , (40) 

where \i is the polarization index of the gluon and a is its color index. This can be approx- 
imated by 

-ig2(l-x B )P B uH a . (41) 

Here = (0, 1, 0) is a vector in the minus direction, so u ■ A = A + . 

What is the net result? The part of the diagram referring to meson B and the interaction 
with the Glauber gluon begins as 

K(k B) q)= (42) 



- iG TT, i « T— (-W)( 2P B ~ - qYt a . 



(k B + q) 2 — m? + ie (P B — k B — q) 2 — m? + ie 
With the replacements described above, this factor becomes 

F£(k B ,q) = - jr - v? { l~ XB l ^7- 2 M<? + ) i^ta. (43) 

(k_b + q) + m z — x B (l — x B )M l v ' 

We can write F as 

F£{k B , q) = -— ^(x B , k B + q) J£(q) , (44) 
x B 

where 

J a ^q) = 2n6(q + ) iguH a . (45) 



and ip is the meson wave function, Eq. 

We can understand this as follows. Let A%(q) represent the lower half of the graph. 
(A depends on other variables, but we suppress this dependence here and think of A as a 
classical gluon field.) Then the factor 

A F a (k B ,q)AZ(q) (46) 



(2tt) 4 ^ 
is approximated by 

/ 4^ ^V>(zB,k B + q) 2n5(q + ) igu^A^q). (47) 
Define Fourier transforms by 
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Kb) = J etq x . 

if)(x B , k B + q) = J d 2 v e'^ B+ ^ r ij(x B , r) . (48) 
Then the factor in Eq. ( PH) can be written as 

/ d 2 r e" ikB ' r il>(x B ,r) / dx~igt a A+{0,x~ , -r) . (49) 

Xb J J-oo 

One can interpret Eq. ( f4"9| ) as follows. The amplitude to find a quark at the annihilation point 
(say x v = 0) with momentum fraction xb while the antiquark carries transverse momentum 
— k B is / d 2 r exp(— ik B ■ r)^(i B ,r). For the present graph, this amplitude is modified by 
multiplication by a line integral of the color potential along a line in the minus-direction 
at the transverse coordinate of the antiquark, — r and the plus-coordinate of the antiquark 
x + ~ 0. The antiquark is not appreciably deflected, but provides a color current that can 
absorb the gluon that was needed to balance the color and maintain meson A as a color 
singlet. The overall minus sign in Eq. (|4"9"D arises because the spectator parton is an antiquark 
instead of a quark. 

spectatorsq.epsf here: just uncomment the macro. 

FIG. 6. Cut graphs for inclusive jet production. Graph a) contributes to diffractive jet produc- 
tion, while graph b) does not. 

We have seen that the contribution to the cross section calculated from the graph of 
Fig. |5] does not take the factorized form of Eq. (|TT|) . To see why the usual factorization 
theorem JTSJ does not apply here, consider the square of Fig. ||, which is represented by 
the cut Feynman diagram of Fig. |6|(a). If we were considering an inclusive jet production 
cross section, then other graphs would also be present and would cancel Fig. [](a) to leading 
order in mjE^. These other graphs do not contribute to the cross section for diffractive jet 
production. For instance, one graph that contributes to inclusive jet production but not to 
diffractive jet production is shown in Fig. |6](b). In this graph, the constituents of hadron A, 
rather than a diffractively scattered version of hadron A, appear in the final state. 



IX. COLOR EXCHANGE WITH ACTIVE PARTON 

We now consider graphs in which the gluon from meson A is absorbed on the active 
quark from meson B before it enters the hard interaction. A typical graph of this sort is 
shown in Fig. [7]. For convenience, we shall refer to the gluon from meson A as the Glauber 
gluon, although the important integration region for its momentum is, strictly speaking, 
not confined to the Glauber region as described in the previous section. We will concentrate 
on the factor associated with the absorption of the Glauber gluon on the active quark line 
from hadron B, 

Gl = 71 Z J TT ~ *9 ( 2kB + ■ ( 5 °) 

[Kb + q) — m + «e 
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active. epsf here: just uncomment the macro. 

FIG. 7. Graph with active parton interaction. 



As in the case of Fig. |5|, the dominant integration region for the transverse components of 
is q 2 ~ m 2 and the dominant integration region for q~ is q~ ~ m 2 / Pj[ . 
The factor ig (2k,B + q) fl t a in ( [50|) multiplies a vector associated with meson A that 
is predominantly in the plus direction. Thus we can replace this factor by ig2xBpBU fl t a 
with xb = k^/Pg and = (0,1,0). We now examine the integration over q + . The 
only propagator that has a significant dependence on q + for small q + is the active quark 
propagator in (|50"D , 

(51) 



(k B + q) 2 — m? + ie 2x b Pb Q + — -AJ + ie ' 

where A| is given in Eq. (|5T|). Recall that A 2 is of order m 2 in the dominant integration 
region. Thus there is a pole in the lower half q + plane near q + = 0. However, as discussed 



in Ref. |[L2| , the integration contour is not pinched between two nearby poles. Thus the 
contour can be deformed away from the pole. We then see that it is a good approximation 
to replace A| by in fl5"i"|) . Then 

Qi « iguH a —^— . (52) 
q + + %e 

activemod.epsf here: just uncomment the macro. 

FIG. 8. Graphs with final state jet interaction. 

Before completing work on the interaction with the active quark entering the hard inter- 
action, we need to consider the two graphs, as shown in Fig. |^, in which the gluon attaches 
to one of the partons emerging from the hard scattering. Here the factor associated with 
the absorption of the gluon is approximately 

s^i^k F - qY t a {kF _ q) ;_ ml+i( , (53) 

where t a is the appropriate color matrix, depending on whether the final state parton in 
question is a quark or a gluon, and where mp is m for a quark or for a gluon. This 
is exactly the vertex factor for the graph in which the Glauber gluon is absorbed on the 
outgoing quark line. It is approximate in the case that the gluon with small momentum q^ 
is absorbed on the outgoing gluon line that carries large momentum kp into the final state. 
As before, the factor ig (2kp — q)^t a can be approximated by ig 2kpU fJ, t a . The propagator 
can be approximated by 



(/cf — q) 2 — m F + ie —2kp ■ q + q 2 + ie 

i 



-2kp ■ q + ie 



i 



-2k F q + - 2kjq- + 2k F ■ q + ie 



(54) 
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This propagator has a pole in the upper half q + plane. The poles in the q + plane from the 
other propagators are all at large q + , q + ~ E?- Thus we can deform the contour into the 
lower half q + plane until |g + | ^> m. Now the components of kp are all of the same order 
of magnitude (about Ep), while the minus and transverse components of q are both small. 
Thus on the deformed contour, the approximation 

(55) 



(hp — q) 2 — mp + ie —2k F q + + ie ' 

is valid. Of course, the integral is independent of the contour deformation, so we can move 
the contour back to the real q + axis. (Similar manipulations may be found in Ref. ||12|| .) 
With these approximations, Qp is 

Q F « -ig U n a —^— . (56) 
q + — ie 

Notice that there is a certain similarity among the effective factor Qj giving the coupling 
of the Glauber gluon to the incoming parton from hadron B and the two factors Qp giving 
the coupling to the outgoing partons. For the moment, let us consider q + ^ 0, so that 
we can ignore the ie. Then the only difference is how the color matrix t a for the Glauber 
gluon is connected to the color matrix M l c 3 d for the hard scattering. Here the color indices 
are i for the final quark, j is the initial quark, c for the final gluon, and d for the initial 
gluon. Attaching to the final quark, we have (t^)u' M*j, where the superscript (3) denotes 
the triplet representation of SU(3). Attaching to the final gluon, we have (t a 8 ^) cc < M l J d . 
Attaching to the initial quark, we have —M l cd {t^)j'j, where we have noted the opposite 
signs between Eq. ( p^) and Eq. (|56|). We can sum the three contributions using the group 
theoretic identity 

(t a 3) h> <l + (*?>)«* M% - M% (^)fj = Mli , (57) 

which merely says that M is invariant under SU(3) rotations. In the color matrix on the right 
hand side of fl57|), the Glauber gluon from hadron A first couples to the active gluon from 
hadron A to make a color octet object, then this color octet couples to the hard scattering 
color matrix. However, these two gluons are in a color singlet state, so coupling the two to 
make a color octet gives zero. (That is, the lower part of the graph is proportional to 5 a d, 
while (t^)d'd is antisymmetric in {a,d}.) 

Thus the sum of the three diagrams cancels, except for the fact that the denominators 
in Eqs. (|5^) and fl56"D have the opposite ze's. In order to take the ze's into account we now 
make the replacement 

+ 27i5(q + ). (58) 



q + + ie q + — ie 

in Qj. The first term has the right ie to participate in the cancellation, and we are left with 
an uncancelled term 

g'j = igun a 2ir5{q + ) . (59) 
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We are grateful to John Collins for suggesting this to us. 

What is the net result? The part of the diagram referring to meson B and the interaction 
with the Glauber gluon begins as 

W'-') = ST? ^«(*+')"'-l|r^«' (60) 

With the replacements described above, this factor becomes 

K"fa,q) = 7T - G J l ' X n 27r5 (<? + ) ( 61 ) 

[Kb + m) — xb{1 — xb)M z v 7 



We can write F' as 



F' a »{k B , q) = — Tp{x B , kfl) JM , (62) 
Xb 



where J and ip are defined as in Eq. 

It is instructive to write this in transverse coordinate space, as we did for the attachment 
to the spectator antiquark. We let A£(q) represent the lower half of the graph. Then the 
factor 



is approximated by 



(27T) 



d^a 1 

' " V0&B,k B ) 2tt 5(q + ) igu^A^q). (64) 



(2tt) 4 x b 

Define Fourier transforms as before, we find that the factor in Eq. fl63|) becomes 

— / d 2 r e~ ikBr $(x b ,t) / dx~igt a A+(0, x~, 0) . (65) 
Xb J J—oo 

This is the same as the factor (|49|) , except for the overall sign and except that the Glauber 
gluon is absorbed at transverse position instead of transverse position — r. 



X. RESULTS 

We have investigated diffractive jet production within the context of a certain simple 
model. We have seen that the cross section has contributions that are lossless in the sense 
that the plus-momentum zP^ transferred from hadron A is transferred without loss to the 
jet system. Thus these terms contain a factor 5(1 — Xj^/z). We now assemble the results 
and rewrite them in a suggestive form. 

We can build up the result using four ingredients. The first ingredient is the line integral 
of the color potential along a lightlike line in the minus direction (the direction of hadron 
B) at a transverse position b, where the hard interaction is at transverse position 0: 
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/+oo 
d y -J2igtcAt(o,y-,b). 
-oo r- 



(66) 



Notice that, in addition to being an operator, A is a matrix that acts on vectors in the 
3-representation of color SU(3). The next ingredient is the color field operator F^(0, 0, 0), 
which destroys a gluon of transverse polarization j and color a at the origin of space-time. 
Combining these operators, we form the amplitude to annihilate the gluon at the origin 
along with one more Glauber gluon at transverse position b, while scattering the meson A 
with momentum transfer (z,t): 



G{(b;t,z) 



1 



4vrX A P+ 



(P;|T{^(b)F^'(0,0,0)}|P A > 



(67) 



(We have indicated a time-ordered product here, but within the context of the model the 
operators commute, so the time-ordering is not relevant.) The next ingredient is the wave 
function to find the quark in meson B carrying a momentum fraction X B and separated 
from the antiquark by a transverse separation r, i/j(X b , r). 

Finally, we have a hard scattering function H^ b (s, Et)jj, where s = XaXbs. Up to a 
normalization, this function is the product of the Born-level scattering amplitude M. for the 
gluon of plus- momentum xaPa to scatter from a (scalar) quark of minus-momentum x b Pb 
times the complex conjugate amplitude Ai> with different color and spin labels. The gluon 
has of polarization k and color b, while the quark has color J. In the complex conjugate 
amplitude, the corresponding indices are j, a, and /. We consider if as a matrix in the 
quark color space and do not write the indices IJ explicitly. The normalization is such that 
the usual hard scattering cross section da for inclusive jet production, as used in Eq. (^), 
is, at the Born level, the color and spin average of H, with the momentum-conserving delta 
functions removed: 



da (a + b — > jets + X) 
dE T dX A dX B 



- Born 



6(x A - X A ) S(x b -X b )WVL \^ { H H} 



With these ingredients, the lossless part of the cross section for diflractive jet production 
within the model described in the preceding sections is 



da dltt (A + B -> A + jets + X) 
dE T dXA dX B dz dt 

2 8 



o 



5(1- X A /z) J dr 



x £ £ Tr^ Gi{-v-t,z)-G{{Q-t,z) 

j,k=l a, 6=1 



mX B ,r)\ 2 
2X B (l-X B ) 



Hit (*, Ei 



x 



G k b (-r;t,z)-G k b (0;t,z) 



(69) 



There are several features of this result that are notable. First is the 6(1 — Xa/z), which 
provides the experimental signature for the lossless contribution to the cross section. In 
contrast, the normal contributions yield values of the variable Xa that are spread over the 
region Xa < z. The second feature is that the factorization that applies to inclusive jet 
production does not apply here. The factor corresponding to hadron A is tied to the factor 
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representing hadron B by a convolution over a transverse position r. There is also a more 
complicated color and spin structure than in the normal contributions. 

We emphasize that the fact that factorization does not apply here does not represent 
a failure of the factorization theorem for the total cross section for A + B — > jets + X. 
The factorization theorem requires an unrestricted sum over the final state X ||12|| . When 
we demand that X contain the a diffractively scattered hadron A, the conditions for the 
theorem are violated. 

Despite its rather complicated structure, the interpretation of Eq. (^) is straightforward. 
In the model, meson B consists of a quark and an antiquark. With probability oc \i/j(Xb, r) | 2 , 
they are separated by a transverse distance r. In order to restore the color of hadron A, 
we must absorb a gluon on either the antiquark (at position — r) or the quark (at position 
0). Since the quark and antiquark have opposite color charges, the absorption amplitude is 
proportional to the difference G{( — r; t, z) — G{(0; t, z). Here G{(b; t, z) is the amplitude to 
absorb a color field quantum at transverse position b when the "active" gluon is annihilated 
at the origin of space-time and hadron A is diffractively scattered. Thus G describes the color 
field associated with the pomeron when one gluon from the pomeron has been annihilated 
at the origin. 

Here we meet an interesting experimental possibility. The b dependence of G J a (h;t, z) 
reflects the transverse structure of the pomeron. It has significant structure on some distance 
scale Rp characteristic of the pomeron. In the present model, 1 / -Rp is of order of the quark 
mass m. Thus G{(— r;t, z) — G{(0;t,z) is small when |r| <C Rp. On the other hand, 
\i[)(Xb, r)| 2 is small when |r| 3> Rb, where Rb is a characteristic size of hadron B. This size 
is also of order 1/m in the model. However, suppose that we generalize the model so that Rb 
can be separately adjusted. Then when Rb ~ -Rp, there will be a substantial contribution 
to the cross section proportional to 5(1 — Xa/z). But when R B <C -Rp, this contribution 
will vanish. 

So far, we have worked only with a simple model. But the model suggests a plausible 
conjecture. First, there can be a sizable contribution to diffractive jet production propor- 
tional to 5 (1 — Xa/z), arising from using one gluon from the pomeron to make the jets and 
absorbing on the partons of hadron B the rest of the color field needed to make hadron 
A back into a color singlet. Second, when the size Rb of hadron B is small compared to 
the transverse size Rp associated with the color field in pomeron exchange, then hadron B 
should act as a color singlet and this contribution should disappear. 

In order to test this conjecture, and probe the transverse structure of the pomeron, 
one needs to use hadrons of adjustable size. This is simple (for theorists). At HERA, 
one manufactures bremsstrahlung photons from the electron beam. The virtuality Q = 
[— PbPbij] 1 ^ 2 of the photon is measured by the deflection of the electron, and can be anything 
from nearly zero to many GeV. The photon can collide with a proton (hadron A) to make jets 
with Ej- ^> Q. The cross section for this process can be (roughly) divided into two parts. In 
one part, the photon acts as a parton and scatters directly with a parton from hadron A to 
make the jets. In the other part, the photon acts as a hadron, made of constituent partons. 
For Q fa 0, this hadron is essentially a p-meson, with a size Rb ~ 1 fm. For Q 1 fm _ , 
the "hadron" consists of a quark-antiquark pair, with wave functions given in Eq. ( p7[ ) for 
transverse polarization and Eq. (|28|) for longitudinal polarization. These wave functions are 
characterized by a size Rb ~ [A B (1 — X B )Q 2 ]^ 1 ^ 2 ■ Since Q 2 and X B are measurable, this 
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size is adjustable. 

Thus the experimental possibility is to select photon-proton scattering events with Q in 
the range from zero up to one or two GeV and look for jets accompanied by a diffractively 
scattered proton. These can be divided into events with a partonic photon (X B 1) and 
those with a hadronic photon (X B < 1). Some of the events should be of the lossless 
(Xa ~ z) type, while some most should be of the normal (Xa < z) type . For the partonic- 
photon events, there should be diffractive jet production, but it should have Xa < z, since 
the photon has no color structure. For the hadronic-photon events, both Xa < z and Xa ~ z 
events should occur, but the fraction of Xa ~ z events should approach zero as Q becomes 
large. The scale of X B (1 — X B )Q 2 at which this happens reflects the transverse size of the 
pomeron. 

In the model we have used, the "lossless" contributions to the cross section are propor- 
tional to 5(1 — Xa/ z) while the "normal" contributions are non-singular functions of Xa/ z. 
We expect that this clean division is not so clean when higher-order contributions are taken 
into account. Presumably there are contributions that are singular as Xa/ z — > 1 but are not 
as singular as a delta-function. Thus one may anticipate that the lowest order 5(1 — Xa/ z) 
appears as events with Xa ~ z rather than Xa = z. 

We must emphasize that the proposal given above is a conjecture based on a simple 
model, not a proven consequence of QCD. It should be a challenge to investigate the structure 
of diffractive hard scattering further and to discover what features of the model survive a 
higher order analysis. 
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